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Measurement of the hyperfine splitting of the 6S1/2 level in rubidium. 

A. Perez Galvan, Y. Zhao, L. A. Orozco^ 
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We present a measurement of the hyperfine spUtting of the 6S1/2 excited level of rubidium us- 
ing two photon absorption spectroscopy in a glass cell. The values we obtain for the magnetic 
dipole constant A are 239.18(03) MHz and 807.66(08) MHz for ^^^Rb and *''Rb, respectively. The 
combination of the magnetic moments of the two isotopes and our measurements show a hyper- 
fine anomaly in this atomic excited state. The observed hyperfine anomaly difference has a value 
of 87585 ~ —0.0036(2) due to the finite distribution of nuclear magnetization, the Bohr-Weisskopf 
effect. 

PACS numbers: 32.10.Fn, 32.30.-r, 21.60.-n 
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I. INTRODUCTION 

Precise measurements of hyperfine splittings of the 
ground and excited states are necessary to explore the 
complete dynamics of the electron cloud-nucleus interac- 
tion in the atom. New experimental approaches such as 
femtosecond frequency combs and small linewidth lasers 
together with laser cooling and trapping reach now in- 
creased accuracy for high precision studies of hyperfine 
structure in excited levels [Ij, t^, ]3-, 4, 5]. This new wave of 
experiments has renewed the interest of theorists in pre- 
dicting accurate electron-nucleus interactions. These cal- 
culations of hyperfine splittings in excited states, where 
electron correlations are less complicated, are more sen- 
sitive to nuclear structure details [Q] . 

Measurements of hyperfine splittings are also of inter- 
est to the atomic parity non-conservation (PNC) com- 
munity. Experiments of atomic PNC rely heavily on cal- 
culations of operator expectation values to extract from 
the experimental data information on the weak interac- 
tion. The accuracy of the calculations is gauged against 
expectations values of atomic properties such as energy 
levels, ionization energy, electric dipole operators of the 
electronic levels, fine and hyperfine splittings. The hyper- 
fine splitting measurements, in particular, represent ideal 
benchmarks for the ah initio calculations of the electronic 
wave function at distances close to the nucleus 0, Q. 
Currently, the PNC ah initio calculations of theoretical 
groups using many body perturbation theory (MBPT) 
have reached a precision better than 1% [ol. [lOl. [l]| . 

We present in this paper the details of the measure- 
ment of the hyperfine splitting of the 6S1/2 excited level 
in ^^Rb and ^^Rb [12|. We perform the experiment in 
a glass cell with rubidium vapor with natural isotopic 
abundances under a controlled environment. We observe 
in our experimental data deviations from the assumed 
point interaction between the valence electron and the 
nucleus i.e. a hyperfine anomaly, and find that a change 
in the distribution of the nuclear magnetization between 
isotopes explains the observation. 

The organization of the paper is as follows: section II 



gives the theoretical background, section III contains the 
methodology and experimental setup as well as the stud- 
ies of possible systematic errors and results. Section IV 
shows the comparison with theory and section V contains 
the conclusions. 



II. THEORETICAL BACKGROUND 

Pauli suggested in 1924 that the "hyperfine" split- 
tings observed in very precise spectroscopic studies of fine 
structure in atoms were due to the interaction between 
the atomic electrons and the magnetic moment of the 
nucleus [1J|. Three years later Back and Goudsmit suc- 
ceeded in analyzing the very small splittings of bismuth 
using the assumption of Pauli and the coupling of angular 
momenta [ij, ll5|. Theoretical predictions of the split- 
tings were moderately successful in explaining the size of 
the experimental findings, mainly due to the complexity 
of the relativistic many-electron system. However, in the 
last twenty years there has been an exceptional output 
of very accurate theoretical results coming from MBPT 
culminating in the extraction of weak interaction cou- 
plings from the atomic PNC measurement in Cs [la, Il3| ■ 
There are currently proposals to perform a complemen- 
tary PNC measurement in francium where calculations 
of atomic properties are reaching the precision of those 
in cesium 11811. 



A. Hyperfine splitting 

Although a complete treatment requires a full rela- 
tivistic theory, estimations of the interplay between the 
nuclear moments and the electromagnetic fields created 
by the electron following classical electrodynamics agree 
with the experimental results and provide physical in- 
sight of the phenomenon (sec for example Refs.|19l.l2C)|). 
We follow this approach in the discussion below. The in- 
terested reader should consult Ref. [2l| for a relativistic 
derivation. 



Two types of nucleus-electron interactions suffice to 
account for the hyperfine splitting in most atoms. The 
largest of the contributions comes from the nuclear mag- 
netic dipole coupling to the magnetic field created by 
the electrons at the nucleus. The second one arises from 
the interaction between the nuclear electric quadrupole 
and the gradient of the electric field generated by the 
electrons at the nucleus. The latter vanishes for spheri- 
cally symmetric charge distributions which correspond to 
electronic angular momentum J equal to or 1/2. The 
hyperfine energy shift Ejjf for these levels is |2g |: 

Ehf = ^{F{F + 1)-I{I+1)-J{J + 1)), (1) 

where F is the total angular momentum, / is the nuclear 
spin and A is the magnetic dipole interaction constant. 
The derivation of A for a hydrogen-like atom by Fermi 
and Segre assumes a point nuclear magnetic dipole [19| 
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where V'(O) is the electronic wave function evaluated at 
the nucleus, jjlb is the Bohr magneton, ^m is the nuclear 
magneton and gj is the nuclear gyromagnetic factor. 

The nuclear magnetic dipole acquires an extra poten- 
tial energy under an external DC magnetic field. For 
small values of the field {gpUBB/EHF <C 1) f is a good 
quantum number and the energy of the system is given 
by 



Ehf{B) = Ehf{^) + gpfJ-smFB, 



(3) 



where gp is the total g-factor, mp is the magnetic quan- 
tum number, B is the magnetic field and E//^(0) is the 
value of the energy at zero magnetic field. In this regime 
of small splittings compared to E/fi?(0), gp is given by: 
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where gj is the electronic g-factor. 



Coulomb interactions with a frozen core Dirac-Hartree- 
Fock wave function of a one valence electron atom. The 
Hamiltonian includes projection operators to positive en- 
ergy states of the Dirac Hamiltonian. Their presence 
gives normalizablc, bound state solutions. The wave 
function contains single and double excitations to all or- 
ders; these correspond to wave functions useful for cal- 
culating energy levels and transition matrix elements. In 
order to calculate accurate hyperfine constants a set of 
triple excitations has to be added. The evaluation of the 
wave function yields coupled equations that are solved 
iteratively for the excitation coefficients which are then 
used to obtain atomic properties. 

The calculations of the hyperfine constants are cor- 
rected for the finite size of the nuclear magnetic moment 
up to zeroth order only due to their small size in the 
lighter alkalies (Na, K, Rb). In cesium and francium 
the correction becomes more important and is included 
to all orders. The calculation ignores isotopic changes 
of the magnetization distribution and it is modeled as 
a uniformly magnetized sphere for all the atoms. The 
magnetization radius used is equal to the charge radius 
and the neutron skin contribution is ignored. 



C. Hyperfine anomalies 

The atomic electron sees the nucleus, most of the time, 
as a structureless entity with a single relevant parameter, 
its charge Z. We should expect the electronic wave func- 
tions of different isotopes, to a very good approximation, 
to be the same. It follows then, using Eq. [2] that 
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where the superindex denotes the atomic number of the 
isotope. 

However, high precision experiments show differences 
or anomalies from this description. It is necessary to 
consider the nucleus as an extended, structured object 
with specific finite magnetization and electric charge dis- 
tributions for each isotope. We can express the anomaly 
by writing the magnetic dipole constant of an extended 



nucleus A^xt as a small correction to A. 
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B. Ab initio calculations 

Hyperfine interactions represent a formidable many 
body problem. A thorough study must approach the 
problem from a relativistic standpoint which further com- 
plicates the interactions in a multi-electron atom. In re- 
cent years relativistic MBPT has shown itself to be a 
powerful and systematic way of extracting, from the high 
quality wave functions that it generates, precise atomic 
properties such as hyperfine splittings [ill . [23 | . 

The full method is outlined in Ref.[2^ and references 
therein. Briefly, the method, applied to alkali atoms, con- 
sists of evaluating a no-pair relativistic Hamiltonian with 



Aext — ApointfR{^ + ('BCRS){^+ ^Bw), 

(5) 

where fji represents the relativistic correction. The last 
two terms in parenthesis modify the hyperfine interaction 
to account for an extended nucleus. The Breit-Crawford- 
Rosenthal-Schawlow (BCRS) correction (li, [H, [H , the 
largest of the two, modifies the electronic wave function 
inside the nucleus as a function of the specific details 
of the nuclear charge distribution. The second one, the 
Bohr-Weisskopf (BW) correction [23, describes the in- 
fiuence on the hyperfine interaction of finite space distri- 
bution of the nuclear magnetization. 
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FIG. 1: Plot of nuclear charge radius of rubidium as a function 
of atomic number. Adapted from Ref. [301 ] . 



Direct extraction of hyperfine anomalies from the ex- 
perimental data requires theoretical knowledge of both 
hyperfine structure constants and magnetic moments. 
However, the anomalies can still be observed from the 
measurements of the magnetic dipole constants in dif- 
ferent isotopes and the values of the g- factors [23, [29^ . 
Deviations from Eq. 3] are expressed in terms of the hy- 
perfine anomaly difference 87^85- 
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with 87(^85 = e*'^ 
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indicates the presence of a hyperfine anomaly. 



1. Breit-Crawford-Rosenthal-Schawlow effect 



nuclear charge. The neutron shell for ^^Rb is closed at 
magic number N — 50 making it impervious to the ad- 
dition and substraction of nuclear matter [30, [3l| ■ The 
substraction of two neutrons to form ®^Rb does not affect 
significantly the electric charge distribution, and the elec- 
tric potential, compared to the one from *^Rb, remains 
the same (see Fig. [T]). 

The expression of ebcrs for the uniformly charged 
sphere and charge on surface models is |2l| : 
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where p is a constant of order unity, p — ^ \<? — {Za)"^, oq 
and a are the Bohr radius and electromagnetic coupling 
constant, respectively, rjv is the nuclear radius, and k 
is related to the electronic angular momentum through 
the equation k = 1 + J{J + 1) - L{L + 1) - S{S + 1). 
Table |I| shows the value of the correction for a uniformly 
distributed charge as well as the nuclear radius of each 
isotope employed in the calculation. 

rjv [fm] Ref. escRS 



^'^Rb 4.2031(18) [30] 0.0090835(34) 
^''Rb 4.1981(17) [30] 0.0090735(36) 

TABLE I: Values of €bcrs and corresponding nuclear radius 
for both rubidium isotopes. 

Rosenfeld and Stroke propose a trapezoidal charge dis- 
tribution to approximate their model. The interested 
reader should consult Ref. 26] for further explanation. 
All three models give relatively large escRS (~1%), how- 
ever, the difference between both isotopes for all models 



is very small: ^bcrs ~ ^ 
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The interaction between an electron and an atomic 
nucleus is precisely described by the Coulomb potential 
when both of them are far away from each other, no 
matter whether the nucleus is a point or an extended 
source. For interactions that require the nucleus and the 
electron to be at close distances, such as the hyperfine 
interaction, an l/r potential is no longer adequate. The 
correction to the electronic wave function due to the mod- 
ified nuclear potential is known as the Breit-Crawford- 
Rosenthal-Schawlow correction. 

Calculations of cbrcs take into consideration how the 
charge is distributed over the nucleus. Rosenthal and 
Breit considered for their calculation the charge to be on 
the surface of the nucleus [2J] . Schawlow and Crawford 
also calculated the change of the wave function except 
they considered the charge to be uniformly distributed 
in the nucleus [25| . Rosenberg and Stroke proposed later 
on a third model to improve the agreement between the- 
ory and experiment: a diffuse nuclear charge distribution 

The neutron and proton shells in rubidium determine 
the deformation as well as the spatial distribution of the 



2. Bohr- Weisskopf effect 

The interplay between nuclear magnetization with the 
magnetic field created by the atomic electrons causes the 
hyperfine splitting in atoms. A natural extension of hy- 
perfine splitting measurements is to compare models of 
nuclear magnetism. 

Nuclear magnetization is described in terms of nuclear 
moments with the biggest contribution coming from the 
nuclear magnetic dipole moment. The assumption of a 
point magnetic dipole gives good agreement between cal- 
culations and experiment, however it does not provide 
the complete picture. Nuclear magnetization has a finite 
volume. The electron wavefunctions of levels with total 
angular momentum J = 1/2 have a bigger overlap with 
the nucleus and are able to experience the subtle changes 
of the spatial distribution of the nuclear magnetization. 
This wave functions need to be modified to correctly ac- 
count for the hyperfine splitting. 

The corrections £bw to the wave functions due to a 
finite magnetization distribution were first computed by 



Bohr and Weisskopf [27| . They assumed a uniformly dis- 
tributed magnetization over the nucleus for their calcu- 
lation with a predicted e^^r — e^^^ that ranges between 
0.11% and 0.29%. The BW correction roughly scales as 
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where the average is taken over the magnetization distri- 
bution, with (r^/r^)^„ = 3/5 for a uniform magnetiza- 
tion. For rubidium this gives a correction of the order of 
0.2%, however it is strongly dependent on spin and or- 
bital states of the nucleons i.e. on the specifics of the nu- 
clear magnetization. Stroke et. al. performed the same 
calculation using a trapezoidal magnetization distribu- 
tion [32| . Their results agree very well with experimental 
information extracted from the ground state; they calcu- 
late a hyperfine anomaly difference of 0.33%. Both of this 
results are independent for the main quantum number of 
the valence electron |19l |. 

The nuclear shell model predicts that the total mag- 
netic dipole moment has contributions from both the pro- 
ton and the neutron shell, each with orbital and spin 
angular momenta [l9| 
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where g'^J^ and gl are the effective nuclear spin and 
nuclear orbital gyromagnetic ratios, respectively, s and I 
are the nuclear spin and nuclear orbital angular momenta 
and the sum is taken over both shells. The ^-factors have 
the values g''Jf='^.l{2) and ^^^^=1.09(2) f33]. 

The magnetic dipole moment in rubidium comes al- 
most entirely from the vector addition of the orbital and 
spin angular momenta of a single valance proton. The 
neutron shell is almost spherical for both isotopes due 
to its closed shell structure and the contribution to the 
angular momentum from the neutron shell is very small. 

The lighter of the two isotopes, ®^Rb, has the valence 
proton in an almost degenerate / orbital with its spin and 
orbital momenta antialigned yielding a value of 7=5/2. 
Adding two more neutrons to the core shifts the energy 
level of the valence proton to the nearby p orbital and 
aligns both momenta giving the known value of 1—3/2. 
Table HIl presents the theoretical prediction of the nuclear 
magnetic moment using Eq. [5] as well as the experimental 
result. It is indeed remarkable that such a simple model 
reproduces closely the experimental results, particularly 
for the closed nuclear shell structure of *^Rb. 

Three main factors make the two stable isotopes of 
rubidium good candidates for observing the BW effect: 
First the different orientation of the nuclear spin of the 
valence proton with respect to the nuclear orbital angu- 
lar momentum. Second, the small relative difference in 
nuclear charge deformation. Third, the change of orbital 
for the valence proton in the two isotopes. 



5Rb 
''Kb 



2.00 
2.64 



1.35298(10) 
2.75131(12) 






TABLE II: Theoretical and experimental values of the nu- 
clear dipole moment for rubidium. 



D. Two-photon spectroscopy 

We use atomic laser spectroscopy to measure the hy- 
perfine splitting in two isotopes of rubidium. Parity re- 
quires a two photon electric dipole transition to reach 
the 6S1/2 state from the 5Si/2 ground state. We increase 
the probability of transition by using the 5Pi/2 level as 
an intermediate step. We develop a theoretical model of 
the two-photon transition that includes the main physi- 
cal aspects of our atomic system (see Fig. ^ based on a 
density matrix formalism. 

Our experimental setup consists of two counter propa- 
gating laser beams going through a glass cell with rubid- 
ium vapor in a small magnetic field. We lock the laser 
at 795 nm on resonance, the middle step to the 5Pi/2 
level, while we scan the 1.324 /itm laser (from here on 
referred to as the 1.3 /Ltm laser) over the 6S1/2 level and 
observe the absorption of the 795 nm laser. The sys- 
tem can be modeled as a three level atom in which the 
on-resonance middle step enhances the excitation to the 
final step and the counter propagating laser beams help 
suppress the Doppler background (see for example Ref. 
[35|). However, numerical simulations show that we have 
to model our system as a five level atom to include its 
main qualitative feature: optical pumping effects increase 
the absorption of the 795 nm laser when the 1.3 fiia laser 
is on resonance. 

Figure [3] shows our simplified atomic model. We have 
neglected the Doppler effects as well as the Zeeman sub- 
levels in order to keep the calculation as simple as possi- 
ble without losing the main qualitatively features of our 
system. Level |1) represents the lower hyperfine state of 
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FIG. 2: Energy levels relevant to our experiment (energy sep- 
arations not drawn to scale). The numbers correspond to 
*^Rb (*^Rb). Straight arrows correspond to the excitation 
lasers, ondulated arrows to decays. 




model reproduces the increase of absorption observed as 
the second excitation goes into resonance. This can be 
explained in the following way. The laser coupling lev- 
els |1) and 1 2), in the absence of the second excitation, 
pumps the atoms to level |5). On steady state there will 
be little absorption due to a very small amount of atoms 
being transferred from |5) to |1). By adding the second 
excitation a new reservoir of "fresh" unexcited atoms ap- 
pears in level |1). Instead of falling to the non-absorbing 
level 1 5), they travel to level |3) and then decay to the 
initial ground state level through level |4). These "fresh" 
atoms will add to the ground state population and in- 
crease the absorption. 



FIG. 3: Energy level diagram of the theoretical model repre- 
senting our system. The 7i,j corresponds to the decay rate 
between levels \i), \j), Oij is the Rabi frequency relating lev- 
els \i) and \j), and ^23 is the detuning from resonance of the 
exctitation laser between levels |2) and |3). 



the 5Si/2 level while |2) is the upper hyperfine state of the 
5Pi/2- The decay rate between the two levels is 72i/27r = 
6 MHz [3a]. We simplify the hyperfine states of the 6S1/2 
level to just one level with decay rate 732/27r — 3.5 MHz 
[37| . The ground and intermediate levels are coupled by 
the Rabi frequency 0112 while the intermediate and the 
excited levels are coupled by 0123 • The remaining two lev- 
els, |4) and |5), represent all other decay channels out of 
the cascade system and the upper hyperfine ground level, 
respectively. The detuning between levels |I) and |2) is 
zero for our experiment, but we let the detuning between 
levels 1 2) and |3) vary as 623- The total population is 
normalized to one. 

We are left with a set of twenty five linear equations 
for the slowly varying elements of the density matrix anm 
after using the rotating wave approximation. These are 

/ XlknCTkk — InkCTnn) + 
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where ujnm — {En — Em)/h is the transition frequency, 
^nm = —^mn IS the laser frequency connecting the levels. 
The damping rate is given by: 

-L nm -- "Z / \ lnk ~r '~fmk)- 
k 

We solve for cri2 leaving the detuning between levels 
|2) and |3) {S23 — ^^23 — ^23) as a free parameter. We plot 
the negative of the imaginary part of ai2 , which is pro- 
portional to the absorption of level |2), as a function of 
623 for several different sets of parameters. Our five level 



Figure |4] shows samples of our simulation. We have 
plotted the absorption of the laser connecting levels |1) 
and 1 2) as a function of the detuning of the second laser. 
Figure [H (a) shows how the absorption increases as the 
second laser goes on resonance while Fig. |4] (b) shows a 
decrease. Both plots have the same model parameters ex- 
cept for the ratio 741/745. This ratio determines whether 
the atom will be lost or return to the cycle. A ratio 
bigger than one pumps atoms preferentially to level |1) 
rather than level |5) which constitutes a fresh reservoir 
of excitable atoms. 
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FIG. 4: Numerical simulation of the absorption of the 795 nm 
laser as a function of the normalized detuning of the 1.3 fim 
laser to level |3) in units of 721. Both plots have the same 
parameters except for the ratio 741/745. (a) Increase of ab- 
sorption with 741/745 — 2. (b) Decrease of absorption with 
741/745 = 1/2. 



III. MEASUREMENT OF THE HYPERFINE 
SPLITTING 

A. Apparatus 

We use a Coherent 899-01 Titanium Sapphire 
(Ti:sapph) laser with a hnewidth of better than 500 kHz 
tuned to the Dl hue at 795 nm for the first step of 
the transition. A Pound-Drever-HaU (PDH) lock to the 
F=1(2)^F=2(3) transition in *'^Rb (^^Rb) in a separate 
glass cell at room temperature stabilizes the linewidth 
and keeps the 795 nm laser on resonance. An HP 8640B 
signal generator acts as the local oscillator for the lock. 
The 795 nm laser remains on resonance for about 40 min- 
utes, much longer than the time it takes to record a single 
experimental trace. 

A grating narrowed diode laser at 1.3 /im with a 
linewidth better than 500 kHz excites the second transi- 
tion. We scan the frequency of the 1.3 /im laser with a 
triangular shaped voltage ramp from a synthesized func- 
tion generator at 4 Hz applied to the piezo control of the 
grating and monitor its frequency with a wavemeter with 
a precision of ±0.001 cm^^. A fiber-coupled semicon- 
ductor amplifier increases the power of the 1.3 /im laser 
before it goes to a large bandwidth («10 GHz) Electro- 
Optic Modulator (EOM). Another HP 8640B modulates 
this EOM. Fig. [5] shows a block diagram of the experi- 
mental setup. 

A thick glass plate splits the 795 nm laser beam into 
two copropagating beams before going to the glass cell. 
The power of each beam is approximately 10 /iW with 
a diameter of 1 mm. We operate in the low intensity 
regime to avoid power broadening, differential AC stark 
shifts and line splitting effects such as the Autler-Townes 
splitting. Both beams are circularly polarized by a A/4 
waveplate. A counter propagating 1.3 /im laser beam 
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FIG. 5; Block diagram of the experiment. Key for figure PD; 
photodiode, P.D.H.: Pound-Drever-HaU, M: magnetic, BS; 
beamsplitter. 
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FIG. 6: Absorption profile of the 6S1/2, F=l and F=2 hyper- 
fine states of *^Rb with sidebands. The big sideband belongs 
to the F=l peak. The small feature on the side of the F=2 
peak corresponds to the second sideband of the F=l peak. 
The glass cell is in a magnetic field of 0.37 G. 



with a power of 4 mW and approximately equal diame- 
ter overlaps one of the 795 nm beams. The lasers overlap 
to a precision of better than 1 mm along 75 cm. The cell 
resides in the center of a 500-turn solenoid that provides 
a magnetic field of 7.4 Gauss/A contained inside a three 
layered magnetic shield to minimize magnetic field fluc- 
tuations 38]. The middle layer has a higher magnetic 
permeability to avoid saturation effects. The solenoid is 
70 cm long and has a diameter of 11.5 cm. We operate 
under a weak magnetic field (_B «1 Gauss) to work in 
the Zeeman linear regime. 

After the glass cell an independent photodiode detects 
each 795 nm beam. The outputs of the detectors go to a 
differential amplifier to reduce common noise. A digital 
oscilloscope records the output signal for different values 
of modulation, polarization and magnetic field and av- 
erages for about three minutes. The order in which the 
absorption profiles are recorded is random. During the 
experimental runs we monitor the current going to the 
solenoid that provides the quantization axis. A thermo- 
couple measures the changes in temperature inside the 
magnetic shield (24 °C) to within one degree. The op- 
tical attenuation for the Dl line at line center is 0.4 for 
^^Rb and about three times less for ®'^Rb. 



B. Method 

We modulate the 1.3 /im laser to add sidebands at an 
appropriate frequency with a modulation depth (ratio of 
sideband amplitude to carrier amplitude) that ranges be- 
tween 1 and 0.1. The sidebands appear in the absorption 
profile at a distance equal to the modulation from the 
main features and work as an in situ scale (see Fig. [6]). 
We measure their separation as a function of the modula- 
tion for values bigger and smaller than half the hyperfine 



splitting. We interpolate to zero separation to obtain 
half the hyperfine splitting (see Fig. [7|) . This technique 
transfers an optical frequency measurement to a much 
easier frequency measurement in the RF range. 

The size of the main peaks depends on the coupling 
strength between transitions; the size of the sidebands 
(as compared to the main peaks) will be determined by 
the strength of the transition and also on the number of 
sidebands simultaneously on or close to resonance. We 
observe under normal experimental conditions that the 
laser sidebands are both close to resonance (the lower 
frequency sideband to the 6S1/2 F=l and the upper one 
to the F=2 transition) when the carrier is around the 
half point of the splitting. The stronger of the transitions 
(F=l) depopulates the 5Pi/2 F=2 level leaving only a 
few atoms to excite with the upper sideband, hence the 
smaller transmission peak for the sideband corresponding 
to F=2. 

We have also observed a much richer atomic behavior 
by changing the laser intensities, polarizations and mag- 
netic field environment of the glass cell. Optical pumping 
moves the atomic population from one level to another 
quite efficiently. This is manifest in how the peaks change 
in magnitude or just switch from an increase of absorp- 
tion to a decrease (see Fig. [8]) just as our very simple 
theoretical model predicts. These effects point out that 
a careful control of the environment is necessary for a 
successful realization of the experiment. 

The transfer of population by specific selection of po- 
larization and magnetic environment can also be used to 
obtain a better experimental signal. There are several 
options to reach the 6S1 /2 level. From the ground hyper- 
fine states we can do AF = 0, ±1 transitions. We find 
that doing the two step excitation in either a a^ : a^ 
or a" : a^ polarization sequence for the 795 nm and 
1.3 /im lasers, respectively, with a AF = 1 for the first 
step increases the amplitude of the signal. By choosing 
this polarization sequence we increase the probability of 
the atom going to the excited state and avoid placing it 



0.015 



^ 0.010 - 



c 
si 

C 

o 



0.005 



Q- 0.000- 

o 
1/1 

< 

- 0.005 - 



0.015 




F=2 



0.0 2.0 4.0 

Frequency [GHz] 
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in a non-absorbing state [4l| . 

We place the rubidium cell in a uniform magnetic field 
collinear with the propagation vectors of both lasers. The 
magnetic field provides a quantization axis as well as a 
tool to probe systematic effects. The hyperfine separa- 
tion is now dependent on the magnetic field strength and 
the alignment with the laser. We measure the hyperfine 
splitting for different values of the magnetic field and 
polarization making sure that the above polarization se- 
quence is always satisfied. We extract the value of the 
splitting at zero magnetic field from a plot of hyperfine 
splitting as a function of magnetic field. 




300 



FIG. 7: Experimental traces that illustrate sideband crossing 
for *^Rb. The larger resonance corresponds to the F = 2 level, 
the smaller one to the -F = 3 level of the 6S1/2 state. The 
dots correspond to the center of the profiles, the point where 
both lines cross corresponds to half the hyperfine separation. 



C. Systematic effects and results 

We study the contributions of several systematic effects 
that can influence the hyperfine separation measurement. 
We analyze the peak shape model for the non-linear fit 
to obtain the separation of the centers of the profiles, 
scan width and scan rate of the 1.3 /im laser, power of 
the 795 nm and 1.3 /im lasers, optical pumping effects, 
magnetic field effects, and temperature. 

A)Peak shape model and non- linear fit. The absorption 
of a Doppler-broadened two level system as a function of 
laser detuning is a Voigt profile. When a multilevel sys- 
tem is considered it is not trivial to write down the func- 
tional form of the absorption of any of the lasers inter- 
acting with the system (see for example Refs. [33, l40|). 
We fit the experimental data to Voigt, Lorentzian and 
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Gaussian functions to find the line centers and compare 
the results for consistency. 

We use the non-linear fit package of ORIGIN™ to 
fit the above mentioned profiles to search for model- 
dependent systematics. ORIGIN™ uses a Levenberg- 
Marquardt algorithm to minimize the residuals given a 
specified error. The program has been used in the past 
by our group to obtain high precision lifetime measure- 
ments [33, |42j . We use the resolution limit of the 8 bit 
analog to digital converter of the scope for these calcula- 
tions which corresponds to 0.5% of the total scale used. 
Lorentzian and Gaussian fits have three variable parame- 
ters to fit for each peak which correspond to the FWHM, 
the line center, the area under the curve plus a single off- 
set for both peaks. Voigt profiles have an extra parame- 
ter which corresponds to the temperature of the sample. 
ORIGIN^" gives the error of each parameter which de- 
pends on the quality of the data. 

Voigt profiles are in very good agreement with the line- 
shape. The fit yields the low temperature limit of the 
Voigt profile i.e. a Lorentzian, and hence is in agreement 
with the linecenter extracted using a Lorentzian profile. 
This is expected since the contribution of the Doppler ef- 
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X^ is 2.13. (c) Normalized residuals of the Gaussian fit, the 
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feet on the resonance lineshape should be minimized by 
the counter propagating laser setup and by an expected 
group velocity selection arising from the the two-step ex- 
citation process ( "two-color hole burning" ) . The 795 nm 
laser will only interact with a small number of group ve- 
locities; these groups will be the only ones that will be 
excited to the 6S1 /2 level by the 1.3 fiia laser. Line cen- 
ters extracted from Gaussian fits agree with results from 
the above mention profiles but decay too fast for fre- 
quencies far away from the centers. We also fit the data 
to a convolution of Lorentzian profiles with a rectangular 
transmission function and an exponential of a Lorentzian 
to search for systematic errors and to understand better 
our residues. 

All peak shape models give line centers consistent 
among themselves. All of them have similar structures in 
the residues within the line width of the resonances (see 
Fig. [H]). We have determined that these features come 
about from the high sensitivity from deviations from a 
perfect fit that a difference of two peak profiles has. In 
other words, by taking the residues we are effectively tak- 
ing the derivative of a peak profile that will be as sensitive 
as sharp the linewidth is. To further verify this we take 
the numerical derivative of the data to search for residual 
structure that might change our measurement. We fit a 
straight line to the data that lies within the linewidth 
and extract when the line crosses zero. The results are 
consistent with the fits. Close analysis of the derivative 
in this region reveals no structure. 

Of the fitted functions Lorentzians yield the smallest 
X^. The fitting error of the line centers for all our data for 
Lorentzian fits range between 15 kHz and 30 kHz. We 
quote the average of all the fitting errors of our data 
in Table IIIII Fig. [5] shows a zoom of the sidebands 
as well as the residues for a Lorentzian and Gaussian 
fits. We extract the line centers with both models; the 
difference in separation for both models is in this case 

\X1 - X2\Lorentzian - \xi - X2\Gaussian = 0.35(68) MHz. 

The reduced x^ of the non- linear Lorentzian fit for all our 
data ranges between 1 and 10 depending on the noise of 
the signal with a x^ average of 2.4 over twenty fits. We 
do not observe changes in the splitting that depend on 
the frequency range fitted around the resonances. 

The relative angle between both copropagating lasers 
induces a systematic shift on the absolute frequency the 
atoms observe due to the appearance of the v ■ ki depen- 
dence on absorption where v is the velocity of the atom 
and ki is the wave vector for either laser. This angle de- 
pendence on the Doppler shift for our system is almost 
the same for both our lasers since the cosine of the an- 
gle between them differs from one by one part in 10'^. 
Furthermore, any residual effect is minimized since we 
measure frequency differences. 

Just like the line shape, analytic expressions for the 
linewidth are difficult to write down. We perform a nu- 
merical simulation of our five level system presented in 
subsection D of the theoretical background in the pres- 
ence of a room temperature velocity distribution. The 
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FIG. 10: Zeeman plot of the hyperfine separation of the 6S1/2 
level of *^Rb with both lasers circularly polarized to better 
than 95% and linear fit. 



resonances show linewidtlis of the order of 30-40 MHz 
wliich are in very good agreement with experimental re- 
sults. 

Distortions of the lincshapc i.e. asymmetries, depend 
on the detuning of the 795 nm laser from resonance. 
These can induce unwanted systematic errors to the mea- 
surement. Numerical simulations show, following Ref [8|, 
that the separation of the hyperfine splitting depends 
negligibly on the detuning from the Dl line. Neverthe- 
less, we look for any asymmetries in the peaks themselves 
and dependence on the direction of scan during experi- 
mental runs. No correlation with these effects is found. 

We interpolate to zero from a plot of distance between 
the center of the sidebands vs. the modulation frequency 
to obtain half the hyperfine separation. The linear re- 
gression coefficients in this plots differ from one at the 
most in 2 parts in 10"*. Typical errors for the crossovers 
amount to about 200 kHz. 

V>)Scan and linearity of 1.3 /jm laser. Non-linearities 
in the piezo driving the feedback grating, hysteresis ef- 
fects as well as a slow thermal drift on the 1.3 /im laser 
can generate undesired systematics in the measurement. 
We look for non-linearities by sending the voltage mon- 
itor of the piezo to a digital scope with an 8-bit reso- 
lution during the experimental runs as well as monitor 
the absorption peaks for asymmetries. Comparison be- 
tween absorption peaks for both types of scan (low to 
high frequency and vice versa) reveals no systematic ef- 
fects. Analysis of the long term drift of the 1.3 /zm laser 
shows a stability of better than 100 kHz over a 5 min. 
period which is longer than the time we need to take a 
single experimental absorption signal. 

G)Power of the 795 nm and 1.3 fim laser. We look for 
systematic dependence on the hyperfine splitting on the 
power of both lasers. We change the power of the 795 nm 
laser from 4 ^W to 40 /iW while keeping the power of the 
1.3 /im laser constant. Low signal to noise ratio and the 
observation of the Autler Townes splitting determine the 



lower and upper boundaries of this interval, respectively. 

The Autler Townes effect predicts a splitting of the 
middle energy level by the on-resonance first step in a 
three level system that is proportional to the square root 
of its intensity [43J . For our typical experimental condi- 
tions the splitting should be less than 4 MHz, very small 
to be resolved with the observed linewidths of the atomic 
resonances (~ 40 MHz). 

The 1.3 /im laser operates very close to its maximum 
power on normal experimental conditions. The power is 
distributed among the sidebands and the main carrier 
depending on the modulation depth. We gradually de- 
crease the power of the 1.3 /xm to half its operating value 
to detect any dependence on the power. We observe no 
correlation. 

D) Optical pumping effects and magnetic field. Opti- 
cal pumping effects are the most delicate of all the sys- 
tematic effects. Both laser beams are carefully polarized 
using appropriate A/4 waveplates and their polarization 
checked with a rotating polarizer in front of a detector 
to better than 95%. The polarization of the lasers as 
well as their alignment with the magnetic field determine 
the relative size of the peaks [mp sublevels) that form 
the resonances of the 6S1/2 hyperfine levels. Compari- 
son of absorption profiles for a set polarization sequence 
for different values of the magnetic field gives qualita- 
tive information of the alignment between the magnetic 
field and the lasers. The positive and negative magnetic 
field orientations in a perfectly symmetric situation, af- 
ter a switch of polarization sequence, should yield the 
same absorption profile. For everyday experimental con- 
ditions (around 1 G) we observe no difference between 
positive and negative magnetic field directions. We see 
broadening of the profiles at magnetic fields twenty times 
larger but no asymmetries. Differences start appearing 
at around 85 Gauss which suggests good alignment be- 
tween the lasers and the magnetic field as well as good 
control on polarization of both lasers. 

The hyperfine separation vs. magnetic field plot pro- 
vides more quantitative information. Eq. (3) states that 
the plot should be linear with no discontinuities as we 
change the value of the magnetic field from positive to 
negative. Our plots show a smooth transition between 
negative and positive values of the magnetic field within 
experimental error. Fig. [TO] shows a sample of our data 
when both lasers are circularly polarized to better than 
95%. 

We monitor the current of the coil generating the mag- 
netic field to detect any fluctuation in the intensity of the 
field. We observe small current noise that manifest into 
fluctuations at most of the order of milligauss. 

E) Temperature. We analyze the position of the ab- 
sorption peaks as a function of temperature of the cell to 
check for related systematic effects such as collision shifts 
for both isotopes. The temperature of the glass cell is in- 
creased from room temperature (23 °C) up to 40 °C using 
a heat tape wrapped around it. While recording data we 
turn off the heating tape to avoid stray magnetic fields 
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Systematic effects 



uI^pIMRz] iyfJp[MRz] 



Optical pumping effects 


< 0.016 


< 0.029 


Power of 795 nm laser 


< 0.020 


< 0.005 


Power of 1.3 /im laser 


< 0.011 


< 0.011 


Atomic density 


< 0.020 


< 0.010 


Non linear fit 


< 0.028 


< 0.023 


B-field fluctuations 


< 0.015 


< 0.025 


Total Systematic 


< 0.047 


< 0.047 



Statistical error 
TOTAL 



0.100 
0.110 



0.160 
0.167 



TABLE IIL Error budget for the hyperfine splitting measure- 
ment 



experimental results in Eq. [6] we obtain a value for the 
hyperfine anomaly difference of §7^85 =-0-0036(2). This is 
less than one percent difference, well beyond the current 
MBPT theoretical calculation accuracy of the hyperfine 
sphttings. 

**^Rb [MHzl **''Rb [MHzl 



vhf 717.54(10) 1615.32(16) 
A 239.18(03) 807.66(08) 



TABLE IV: Hyperfine splittings {uhf) and magnetic dipole 
constants for the 6S1/2 level. 



generated by the current going through it. The temper- 
ature of the glass cell is monitored with a thermocouple 
inside the magnetic shield with an accuracy of one degree 
celcius. No dependency on temperature is found. 

We have concluded after close analysis of these studies 
that, to the accuracy of our measurement, Gaussianly 
distributed statistical fluctuations dominate our exper- 
iment (see Table IIIip . The statistical error in the hy- 
perfine splitting, as stated by the standard error of the 
mean, is 110 kHz for ^SRb and 167 kHz for ^^Rb. 

Figure [TT] shows the values of the magnetic dipole con- 
stant for *^Rb for all experimental runs of this work. The 
final result of each run is determined by an interpolation 
to zero magnetic field as a function of the current in the 
solenoid. 

Table IIVI contains the measurements of the hyperfine 
splitting of the 6S1/2 level as well as the corresponding 
values of the magnetic dipole constants for both isotopes. 

The precision of our data allows us to observe a hyper- 
fine anomaly. We use the values of Ref. [44] for the ratio 
^85/^87 ^ 0.295055(25). This is consistent with the ex- 
perimental values of Ref. [3J| . Using this value and our 
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FIG. 11: Results of different runs of the magnetic dipole con- 
stants of the 6S1/2 state of *^Rb. The dashed line corresponds 
to the mean, the solid lines to the 1-a error. 



IV. COMPARISON WITH THEORY 

We compare in Figs. [T^and[T51the results from this ex- 
periment with the previous experimental results of Gupta 
et al. [43| and the theoretical predictions of Safronova et 
al. 9] . The hyperfine anomalies are still not within reach 
of MBPT so the value of *^Rb comes from considering 
no hyperfine anomaly. 





SD [MHz] SDpT [MHz] 


Experiment [MHz] 


5Si/2 


642.6 


1011.1 


1011.910813(2) [46] 


5Pl/2 


69.8 


120.4 


120.499 (10) [1] 


6S1/2 


171.6 


238.2 


239.18(3) (this work) 


6P1/2 


24.55 


39.02 


39.11(3) [3, J 


7Si/2 


70.3 


94.3 


94.658(19) [5] 



TABLE V: Single Double (SD) and partial triple (SDpT) ex- 
citation calculated from ab intio MBPT in Ref. [3] and ex- 
periment magnetic dipole constants for the first J=l/2 levels 
in »5Rb. 

Table |V] shows the values of the magnetic dipole con- 
stants using relativistic MBPT Q with single double 
(SD) and single double partial triple (SDpT) wave func- 
tions and values extracted from measurements of the hy- 
perfine splitting in other electronic states currently in the 
literature for J=l/2 [E i, i, i, El] . We have not been 
able to find in the literature values for higher levels with 
adequate precision to include them in the figure. The 
agreement of the theory with the experiment, for J=l/2 
levels, is well within the 1% level. The SDpT relativis- 
tic wave functions do indeed improve the accuracy of the 
calculations of the single double wave functions. 



V. CONCLUSIONS 

We have measured the hyperfine splittings of the 6S1/2 
level of *^Rb and *^Rb to a precision of 103 ppm and 153 
ppm, respectively. Our measurement is consistent with 
and decreases the uncertainty of the past measurements 
[45| by a factor of 63 for ^^Rb and by a factor of 30 for 

85Rb. 
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FIG. 12: Comparison between experimental and theoretical 
results of the magnetic dipole constant of the 6S1/2 state of 
*®Rb. The points labeled "a" and "b" correspond to our work 
and the work of Ref. |45] , respectively. Point "c" corresponds 
to the theoretical prediction of [9|. 
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The dashed line corresponds to the theoretical prediction for 
a diffuse magnetization distribution. See text for references. 



We are able to extract the hyperfine anomaly with our 
experimental data and show that precision measurements 
of the hyperfine structure in atomic states with different 
radial distributions can give information on the nuclear 
magnetization distribution. The hyperfine anomaly dif- 
ference we extract for the 6S1/2 is sT'^ss = —0.0036(2). 
The difference in the anomalies is indeed a factor of thirty 
larger than the expected BCRS contribution and it comes 
from the BW effect. Fig. [TJ] shows that the anomaly 
measured with the nSi/2 levels is the same independent 
of the principal quantum number as well as the smaller 
deviation from thepoint interaction, if any, for the nPj^ /2 
levels 0, i, i, i, Ea Si- Table |Vll shows the hyperfine 
anomaly differences for the first J ~ 1/2 levels. 
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5S1/2 


-0.00356(8) [46] 


5Pl/2 


-0.0001(1) [1] 


6S1/2 


-0.0036(2) (this work) 


6P1/2 


0.0000(8) [3, J 


7Si/2 


-0.0032(2) [5] 



TABLE VI: Hyperfine anomaly differences gr^ss for the first 
J=l/2 levels in rubidium. 



The plot and the table confirm that it is the electronic 
angular momentum and not the principal quantum num- 
ber that determines the Bohr-Weisskopf effect and the 
bigger anomaly of the nSi/2 levels due to their larger elec- 
tronic density near the nucleus [23l ■ These new measure- 
ments invite new calculations of atomic properties and 
constrain nuclear calculations. As the nuclear charge and 
magnetization distribution are better understood they 
will further test and refine the calculations which are 
of crucial importance for parity non-conservation exper- 
iments. 



FIG. 13: Comparison between experimental and theoretical 
results of the magnetic dipole constant of the 6S1/2 state of 
^'^Rb. The points labeled "a" and "b" correspond to our work 
and the work of [4^, respectively. Point "d" corresponds to 
the theoretical prediction of Ref. [0]. The value for *^Rb is 
obtained from the theoretical prediction for *^Rb by consid- 
ering the hyperfine anomaly as nonexistent. 
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